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(Dated: February 7, 2008) 

Abstract 

The K + — > 7r°7r°e + ^ (ifg4 ) decay has been measured with stopped positive kaons for a data sam- 
ple of 216 events. A comparison of the observed spectra with a Monte Carlo simulation determined 
the K®® form factor. The results are compatible with the K + — > ir + ir~e + v data, as estimated from 
the AI = 1/2 rule. We also established that the K®® channel can be used to determine the tt-tt 
scattering lengths. 

PACS numbers: 13.20.Eb, 13.75.Lb 
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K — > Tnrev (K e i) decays are the subject of continuing interest for a variety of reasons Q, 
The decay form factors serve as constraints on the parameters of the Lagrangians 
in the framework of Chiral Perturbation Theory (ChPT). More importantly, they are clean 
sources of tt-tt pairs at low energy, from which one can deduce the scattering lengths of 
significance for models of hadron dynamics. 

Recently, the BNL-E865 group reported on the precise measurement of the K + — > 
7T + 7T~e + v (K^) decay, yielding the decay form factors as well as the S-wave, isospin zero 
7r-7r scattering length |4J. In addition to the K^f decay, there are two other K ei decays, 
K + — > 7r°7r°e + z/ (Kei) and Kl — > 7r~7r°e + z/ (K~®) Since these K ei decays are theo- 
retically related to each other by simple isospin arguments comprehensive experimental 
studies for all channels provide deep insight into the understanding of K ei physics. Among 
the three K e ± decays, the K®2 channel is the simplest because the decay kinematics can be 
described by only one form factor in view of the identity of the two 7r°s in the final state. 
Thus far, only 37 events have been observed using an in-flight-if 4 " technique [l,!^. A 
more accurate measurement of the K®2 channel was desirable as an additional independent 
check on the tt-tt scattering length and the AI = 1/2 rule. 

In this paper, we present a new measurement of K®2 decay. The experiment used a 
stopped K + beam in conjunction with a 12-sector iron-core superconducting spectrometer. 
We have been able to make a first attempt to deduce the q 2 dependent terms of the form 
factor and determine the tt-tt scattering length in a nearly model independent manner. 

The experiment was performed at the KEK 12 GeV proton synchrotron. The experimen- 
tal apparatus was based on the E246 experiment, a search for the T-violation in K + — > tt°li + v 

nn 

decay (K^) [allOj. Besides the T-violation search, spectroscopic studies for various decay 
.nodes W bee. s— y pe— Q. After th e E 246 d a t a eolleefto. was completed 
in 2000, the E470 experiment to measure direct photon emission in the K + — > 7r + 7r°7 (K^-y) 
decay was carried out fljjj], and the events were simultaneously recorded. 

A separated 660 MeV/c K + beam was stopped in an active target system located at the 
center of the spectrometer. The iT™ events were identified by analyzing the e + momentum 
with the spectrometer and detecting the four photons in the CsI(Tl) calorimeter. Charged 
particles from the target were momentum-analyzed using multiwire proportional chambers, 
as well as an array of ring scintillators. The e + s were separated from /i + s by determining 
the mass squared (M| OF ) of the charged particles from the measured time-of-flight between 
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T0F1 and TOF2 counters. Counter TOF1 surrounded the active target and counter TOF2 
was located at the exit of the spectrometer. The photon energy and hit position were 
obtained, respectively, by summing the energy deposits and taking the energy-weighted 
centroid of the crystals sharing the shower. 

K®® decays at rest were extracted by the following procedure. The K + decay time, 
defined as the e + signal in the TOF1 counter, was required to be more than 1.4 ns later 
than the K + arrival time as measured by the Cerenkov counter. This reduced the fraction 
of K + decay in-flight contamination to the level of 10~ 3 of the stopped K + s. Events with 
e + scattering on the magnet pole faces were eliminated by requiring the hit position in the 
ring counters to be consistent with the charged particle track. The selection of positrons 
required —5000 < M| OF < 5000 MeV 2 /c 4 . Events with four photon clusters coming from 
two 7r°s were selected, while events with other photon cluster number were rejected. Photon 
conversion backgrounds in the active target system could be removed by selecting only events 
in which there was one charged particle hit in TOF1 as well as in the ring counters. 

Since there are three possible combinations to form two tt°s from the four photons, the 
K®1 events were reconstructed by introducing a quantity Q 2 , 

Q 2 = (M n o - M 1 f/a 2 Ml + {M A - M 2 ) 2 /a 2 M2 

+ (cos0X" - cos^to - 5) 2 /al (1) 

where M n o is the invariant mass of the selected pair (the first 7r° has higher energy) and 
is the opening angle between the 7r°s. The superscripts MEAS and CALC stand for the 
measured angle and the angle calculated from the measured photon energies, respectively. 
The photon pair with the minimum, Q 2 ^, among the three combinations was adopted as the 
correct pairing. The a and offset values in each term are <jm x = 12.74 MeV/c 2 , o"m 2 — 15.42 
MeV/c 2 , M l = 124.5 MeV/c 2 , M 2 = 113.5 MeV/c 2 , 5 = -0.019, and a 5 = 0.336. The 
choice of Eq. (0) and these parameters were determined to maximize the probability for the 
correct pairing using a GEANT based Monte Carlo code. Because of the shower leakage 
from the calorimeter holes, the a and offset values differ from the ideal ones. The correct 
pairing probability was estimated to be 96% from the simulation. 

Since backgrounds due to K + — > 7r°e + z/ (K e ^) and K + — > 7r°e + z/7 (K e ^) decays with 
accidental photons, photons split into multiple clusters, and bremsstrahlung photons do not 
satisfy the kinematics, the cut of < 3 essentially removed these backgrounds. The 
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fraction of the accidental background in the calorimeter was estimated to be 1.0% by chang- 
ing the CsI(Tl) TDC gate widths. Using the Monte Carlo simulation, background fractions 
in the K ei sample due to the bremsstrahlung photons and photons split into multiple clus- 
ters were obtained to be 0.02% and 0.06%, respectively. The total level of the backgrounds 
was thus derived to be 1.1% which was dominated by K e3y with accidental photons. The 
opening angle distributions between e + and four 7s (0 e + 7 ) and between each photon (6*77) 
for the selected K e ^ events were compared to the simulation to confirm the correctness of 
the simulation conditions, as shown in Fig. ^a)(b). No deviation of 9 e + 1 and # 77 from the 
simulation at small angles was observed, supporting the correct estimation of these back- 
ground fractions. After these selection conditions, the number of good K®® events was found 
to be 216. The e + momentum and 7r° invariant mass spectra of the selected K®® events are 
shown in Fig. fTTc)(d). 
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In general, the K e4 kinematics can be written as a set of five configuration variables |1|, |3J : 
(I) s w = M^, the effective mass squared of the dipion system, (II) si = M^ u , the effective 
mass squared of the dilepton system, (III) 9 n , the angle made by the tt° in the dipion center 
of mass with respect to the kaon direction, (IV) 61, the angle made by the e + in the dilepton 
center of mass with respect to the kaon direction, and (V) <p, the angle between the decay 
planes formed by the dipion and dilepton systems. Because of the identity of the two 7r°s, 
9 n and <fi can be defined in the region of < 9^ < tt/2 and < <p < tt /2. The most general 
form of the K e 4 matrix element in terms of the hadronic vector and axial vector current 
contributions and is given by 

M = ^V us u(Pu)l,(l - IsWPeW - A»), (2) 
A" = FP» + GQ^ + RL^, (3) 
V» = He^ upa L u P p Q a , (4) 

where P = p\ + p 2 , Q = Pi — P2, L = p e + p v , and pi, p 2 , p e , Pu are the four- momenta of 7Ti, 
7T2, e + , v in units of uik, respectively. The form factors F, G, R, and H are dimensionless 
functions of s n . The R term enters the decay distribution multiplied by ml and is therefore 
negligible. In general, S and P waves in the dipion system were considered to contribute 
to the K e 4 decays 0. However, odd angular momenta n-n states cannot be present in the 
^4 channel due to symmetry arguments. Therefore, only a contribution from the S"-wave 
component is possible and the decay kinematics are described only by the F form factor 



Here, a parameterization 



F = fo + f'q 2 + fV, 



(5) 



is employed, where q 2 is the momentum transfer squared q 2 = s n / (4m^.) — 1 to the lepton 
current. The observed spectra (s n , cos^o^o, 6^, 4>) which are sensitive to the K®2 form 
factor are shown in Fig. 

The f'/fo and f"/ fo values could be extracted by fitting the observed spectra to the Monte 
Carlo simulation. The simulation data were generated according to Pais and Treiman |2] and 
analyzed in the same manner as the experimental sample. The simulation spectra were then 
obtained as a function of the form factor by re- weighting the reconstructed simulation events 
using the generated Monte Carlo kinematic variables. A program based on MINUIT mini- 
mized x 2 an d estimated the error of the parameters. The reproducibility of the experimental 
conditions in the simulation was carefully checked using K n2 and K n3 decays. 

In the case of the decay, the tt-tt scattering length can be determined from the 

phase-shift difference between the S and P waves |4j. However, a similar procedure would 
not be of use for the K®2 decay. Instead, one can assume that the s n dependence of F is 
described by an 5-wave phase-shift parameter 5® 



where @ — (1 — 4m^/s 7r ) 1 / 2 and the dependence of the phase shift is given by the Chew- 
Mandelstam effective-range formula [l4| : 



To determine the a[j value, this F form factor was used to generate the K®2 events in the 
simulation. Since the f'/fo and f" / fo parameters are numerically determined from a° by 
fitting Eq. JBJ) to Eq. (J3J), it is possible to check the consistency of the theories by comparing 
the results of the (f'/fo, f'/fo) and a° analyses. 

The f'/fo and f'/fo values were determined from the density distribution of the two 
observables, s n and 6 n o n o. The experimental distribution p(s n ,8 7T o n o) was fitted to a sim- 
ulation spectrum with f'/fo and f'/fo being free parameters. They were obtained to be 
/'//„ = -0.45+°;^ and f'/fo = 0.32l£g. The reduced \ 2 was 0.91. Fig. H shows the 



F oc (l/f3)sm[5^ S7r )}e<^\ 



(6) 



cot5 ° = l/((3a° ) + (2/7r)ln[Vi;/(2m 7r )(l + (5)}. 



(7) 
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X 2 contour plot in the (f'/fo, f"/fo) space. The histograms in Fig. |21 shows the sim- 
ulation spectra with the best fit parameters, which are in good agreement with the ex- 
perimental data. Using these form factors for the acceptance calculation, we obtained 
T(K^) = [(3.56 ± 0.26(stat.)tlH(syst.)} x 10 3 s" 1 by normalizing it to the theoretical 
value of internal bremsstrahlung in the decay [15]. The large systematic error from 

the uncertainty of the form factors (i.e., uncertainty of the e + momentum distribution) was 
due to the fact that our experiment covered a very narrow region of phase space and the 
error of the detector acceptance became large. As shown in Fig. El the experimental 6 n and 
distributions were well accounted for by the simulations without any D-wave component. 
Thus, the data are compatible with an S"-wave only description. The reduced \ 2 was 0-92 
for 9k and 1.19 for <fi. This feature also supports the previous K^ 4 ~ measurement in which 
the D-wave contribution was assumed to be negligible 

Then, using Eq. © as the F form factor, the a[j value which gave the best fit to the 
experimental p(s 7r , cosO^^o) distribution was determined to be a® = 0.45 ± 0.43 in units of 
m n . The reduced x 2 was 0.91. While the large error precludes a meaningful comparison with 
results from other channels, a well designed experiment at the future high intensity facilities 
such as J-PARC will be able to make a significant contribution to the determination of the 
a[j parameter. The f'/fo and f'/fo parameters were determined from this a[j value to be 
/'//„ = -0.301°;! and /"//„ = 0.12±g;i|, which are in agreement with the (/'//„, /"//„) 
fitting. 

The AI =1/2 rule predicts that the q 2 dependence of K®2 and form factors are 

identical. As seen in Fig. 01 our K®2 result overlaps with the K^T form factor reported by 

the BNL-E865 group within two standard deviations |4|. The f'/fo values were obtained 

to be -0.19lg;J? and -O.lll&l? for /"//„ = and -0.10 (K+T result), respectively, which 

also overlaps with the K^ 4 ~ result within the error |4|. The AI =1/2 rule also predicts 

a simple relation for the decay widths r(K° 4 °) = 1/2T{K+-) - l/4r(fQ°) [6]. The world 

averages for the K+ 4 ~ and K~ 4 ° channels [lj give T(K^) = (1.40 ± 0.04) x 10 3 s" 1 . Using 

the form factor j4| for the acceptance calculation of the present experiment to reduce 

the systematic error, we deduced T(K^) = [1.85 ± 0.13(stat) ± 0.2A(syst.)) x 10 3 s _1 . The 

average of the previous K®2 experiments, under the assumption of constant form factor, is 

I I 

quoted as T(K^) = (1.70 ± 0.32) x 10 3 s -1 [16]. While these numbers are all consistent 
with each other, the large systematic error of the T(K^) value determined in the present 
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experiment precludes a meaningful comparison with other results. Further measurements 
will be essential before we can draw firm conclusions on the A J = 3/2 contributions. 

In conclusion, we have performed a new measurement of one of the K e4 channels K + — > 
7r°7r°e + z/ (K^) using stopped positive kaons. The data sample of 216 events is, while a factor 
of six improvement over the world data set, still very limited. The detector response and 
acceptance functions were evaluated by the Monte Carlo simulation. The backgrounds due 
to K e3 and K e3y decays were estimated to be 1.1%. The F form factor of the decay was 
determined by fitting the observed p(s n , O^o) distribution to the simulation. A first attempt 
to deduce the S'-wave n-n scattering length from the p(s n , O^o^o) distribution was made. The 
good reproducibility of the experimental 8 n and <fi spectra by the simulation indicates that 
the data are compatible with an S'-wave only description. These results are consistent with 
the previous K + — > 7r + 7r~e + z/ measurement within errors. By this measurement, despite the 
limited statistics, we established the physics potential of the channel. 
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FIG. 1: -K"^ spectra: (a) cos# e + 7 , (b) cos6> 77 , (c) e + momentum, (d) M n o for the experimental data 
(dots) and a Monte Carlo simulation of K®® (histogram). Four and six combinations for cos# e + 7 
and cos# 77 , respectively, were accumulated in the same figures. Both M n o and M n o values were 
also accumulated in (d). 
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FIG. 2: K®2 spectra: (a) s n , (b) cos^o^o, (c) 6 W , (d) (/>, for the experimental data (dots) and the 
Monte Carlo simulation of K®® with the best fitted f'/fo and /"//o parameters (histograms). 
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FIG. 3: x 2 contour plot in the (f'/fo,f"/fo) space. Contour levels of 1-, 2-, and 3a constraints 
are drawn. The experimental result by the BNL-E865 group is also plotted. 
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